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A possible effect of direct CP violation in D — > Kgiv'^Tv' decay on the 7 measurement from 
tt'^-r Dalitz plot analysis is considered. Systematic uncertainty of 7 coming 
from the current limits on direct CP violation in D — > K%-k'^-k~ is estimated, and a modified model- 
independent procedure of — >■ DK^, D — >■ K^tv'^tt^ Dalitz plot analysis is proposed that gives 
an unbiased 7 measurement even in presence of direct CP violation in charm decays. The technique 
is applicable to other three-body D decays such as D'^ — >■ KgK'^ K~ , — >■ n^n'n", etc. 
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I. INTRODUCTION 

The mechanism of CP violation in particle physics is of 
primary importance because of its impact on cosmologi- 
cal baryogcncsis and possible antimatter existence in the 
universe. In the quark sector, CP violation is studied by 
measuring the elements of Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix with the convenient repre- 

sentation given by the Unitarity Triangle (UT), the an- 
gles and sides of which are parameters observable in var- 
ious decays of B mesons. Precision measurement of the 
UT angle 7 (also denoted as ^3) is an essential ingredi- 
ent in searches for New Physics phenomena in B decays. 
The value of 7 acts as one of the Standard Model ref- 
erence points against which other measurements of the 
UT parameters are compared. The irreducible theoret- 
ical uncertainty in the extraction of the angle 7 from 
B^ — > DK^ decays is due to electroweak corrections 
and is extremely small, of the order of 10~^ How- 
ever, the experimental determination of 7 value remains 
a challenge owing to the low probabilities of the decays 
involved. 

The types of measurements that dominate 7 sensitivity 
are based on B^ — > DK^ decays where the neutral D 
meson decays into a CP eigenstate (commonly referred 
to as the GLW method suppressed Kn state (ADS 

method or self-conjugate three-body final state such 
as K^TT+n- (GGSZ or Dalitz plot method 0,i|). None 
of these methods are systematically limited at the current 
level of precision. However, obtaining a degree-level pre- 
cision on 7 will require some subtle effects to be taken 
into account. The effect of charm mixing has already 
been considered by several authors [9l-[ll| and was found 
to be negligible in most cases. Recent evidence of direct 
CP violation in singly Cabibbo-suppressed two-body D 
decays reported by LHCb has triggered discussions of 



the possible effect of CP violation in charm on 7 measure- 
ments using the GLW technique [Tsl - flH . Although the 
evidence of large CP violation in D hh decays is not 
supported by the updated LHCb measurements [H, |T3| , 
this effect can play its role in precision measurements of 
7- 

In this paper, we consider a possible effect of direct 
CP violation in charm on the measurement of 7 using 
the Dahtz plot analysis of B^ DK^, D K^tt+tt- 
decays. The decay D — )■ Kgir^TT" is dominated by 
Cabibbo-favored transitions, and thus direct CP asym- 
metry coming from the Standard Model effects is ex- 
pected to be very small. However, future precision mea- 
surements of 7 can reach the point where this contribu- 
tion will become significant. On the other hand, if dis- 
agreement in the UT parameters due to New Physics will 
be found in future measurements, the method that can 
distinguish whether the New Physics contribution enters 
charm or B decays will be essential. In addition, "effec- 
tive" CP violation in D /C^tt+tt" decay of the order 
of 10~^ should arise from the CP violation in the neutral 
kaon system if this effect it not explicitly accounted for. 



The goal of this paper is twofold. First, we estimate 
the systematic uncertainty on 7 coming from the current 
limits on direct CP violation in the D -> A'^tt+tt" decay. 
Second, we show that in the model-independent analysis 
using quantum-correlated DD data at charm threshold 
it is possible to account for the CP violation in charm 
and obtain an unbiased measurement of 7 without sig- 
nificantly sacrificing the statistical precision. 

Although the decay D AT^tt+tt^ is used throughout 
this paper, the same approach can be applied to other 
three-body final states of the neutral D decay, such as 



K-, D'> 



7r°, etc. 
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II. FORMALISM OF MODEL-INDEPENDENT 7 
MEASUREMENT WITH CP VIOLATION IN 

D -> K^n+n- 



The procedure to extract 7 from — > DK^, D — s- 
TT^ in a model-independent way employed in the 
current analyses [T8l - [20| assumes CP conservation in D 
decays. The technique uses binned Dalitz plot distri- 
butions, and in order to utilize the assumption of CP 
conservation the bins are chosen symmetrically to the ex- 
change of Dalitz plot variables of the D KgTT^TT~ de- 
cay (invariant masses squared of KgTr~^ and KgTr~ combi- 
nations): rn^o^+ "^/s-Ott^" bins are denoted with 
the index i which runs from —J\f to J\f excluding zero; 
symmetric bins have the same |i| and the flip of the sign 
i —i corresponds to the reflection w^o^+ ^ 'm^j^o^-- 
Current analyses [T8l - [20| are performed with TV = 8: this 
number of bins, together with a special choice of the 
shape of the bins over the Dalitz plot, provides a sta- 
tistical precision for the 7 measurement that approaches 
the precision of the unbinned model-dependent technique 
with the limited DK^ data available today [2l| . 

The procedure of model-independent Dalitz plot anal- 
ysis is described in detail in Refs. 0, [2l|; here we give 
only the final equations. The analyses use four cate- 
gories of events involving D K'^tt^tt^: flavor-tagged 
D — ivTg 7r+7r~, neutral D mesons tagged in CP eigen- 
state from ?/'(3770) — >■ DD process [Dcp), correlated 
pairs of neutral D mesons where both D are recon- 
structed in the /C^tt+tt" state, and B^ DK^ decays 



with D ~» K'^TT+TT- ' . The number of events in bin i of 
the flavor-tagged D decay is different for and 13°; it 
is denoted as Ki and Ki, respectively. However, in the 
CP-conserving case and with the symmetric binning de- 
scribed above Ki ~ K^i, so Kj, are not independent. 
The numbers of events in bins are then related as 



= hcp [K^ + K-^ + 2^K^K-,C^ 
for D decays into a CP eigenstate. 



2^K~K~K~K~{C^CJ + S^SJ) 



(1) 



(2) 



for correlated DD pairs both decaying into K'qTt^ti^ 
(here one has to deal with two correlated Dalitz plots 
and thus the number of events is described with two in- 
dices i and j), and 



NT = h 



2y/KJ{~{x±C^±y±S^) 



for D 
rs cos((5_B ±7), y± 



from B^ — !> DK^. Here x± 



(3) 



tb cos{6b ±7); 



+ yi. The 



free parameters x± and y± hold the information about 
the phase 7 and hadronic parameters in B^ — >■ DK^ 
decay: the amplitude ratio and the strong phase dif- 
ference Sb- The other free parameters are the normal- 
ization factors hcp, /icorr and hB±, and phase terms Ci 



and Si. The terms Ci and Si describe the average sine 
and cosine of the strong phase difference between D'~' and 
D'^ amplitudes over the bin i. In the case of CP conser- 
vation they satisfy Ci — C-i, Si = —S-i. Thus, these 
parameters are independent only for i > 0. The system 
of equations (dl, ©, and is overconstrained and can 
be solved with the maximum likelihood fit to obtain x± 
and y± and, thus, the value of 7. 

Now we turn to the case when CP is not conserved 
in the D — >■ ifg7r+7r~ decay. The relations between 
symmetric bins of the Dalitz plot do not hold any- 
more. We still use the notation for bin number i = 
-A/",... - 1,1,... A/", but now Ki ^ K_„ C, ^ C_^, 
and Si ^ —S^i in general. In principle, now the binning 
is not required to be symmetric, although in our studies 
we keep the same binning as in the CP-conserving case 
to allow for a direct comparison of the two approaches. 

The equations relating the numbers of events in bins 
of the D -I- KgTT^TT^ Dalitz plots are: 



M, 



T-CP 



K, + K, + 2\ K,K,C. 



Mij = /icon- [KiKj + KiKj- 



2JK^K^KjK,{C^Cj + S^Sj 



(4) 



(5) 



and 



K, + rlK, + 2^K,K,{x+C^ + y+S^) 



K, 



K, + 2^jK,K,{x^C^-y-S, 



(6) 



Note that the number of phase terms Ci , Si is doubled 
compared to the CP-conserving case since their values for 
i < are now independent. The numbers of flavor-tagged 
events Ki and Ki also have to be obtained independently, 
but since the available samples of flavor-tagged D decays 
are large, this should not limit the accuracy of the mea- 
surement. There are 4Af +8 free parameters for 4J\f'^+6J\f 
equations (P), ©, and ([3]), and thus the system of equa- 
tions still remains solvable. As a result, arbitrarily large 
CP violation in D decays does not lead to a bias in the 
measurement of the x,y parameters, and, hence, of the 
value of 7 when using this technique. 

The decays of a neutral D m a CP eigenstate into 
KgTT^TT^ are obtained from the process -0(3770) ^> DD, 
where the other (tagging) D meson is reconstructed in 
the CP eigenstate of the opposite parity. Therefore, 
if CP is violated in the decay of the tagging D, Eq. 
(jl]) would not be valid. This effect is expected to be 
larger for CP-even tags using Cabibbo-suppressed decays 
(Z?*^ — > X+X^, TT+TT^) than for CP-odd tags which are 
mostly Cabibbo-favored (such as D° Kgw'^). How- 
ever, even without CP-tagged D decays, the remain- 
ing equations © and which do not include D de- 
cays other than A'^tt+tt", can still be used to perform 
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the measurement. The only limitation arising from this, 
apart from the reduced sensitivity to the parameter C.j, 
is that the system ([5]) , ([6|) has one more additional ambi- 
guity compared to the case when Dcp decays are used: 
it is insensitive to the simultaneous change of sign of all 
Ci. The choice between the two solutions can, though, 
be made using Eq. ^ with the good assumption that 
CP violation in the tagging D decay is small. 

The technique described above can be applied not only 
to D ^ KgTT~^Tr~ decays with CP violation, but also 
to other non-self-conjugate final states, such as I? — ?> 
K%K-'K+, D K+n-n", etc. 



III. BIAS OF 7 MEASUREMENT DUE TO CP 
VIOLATION IN A'g7r+7r DECAY 
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In this section, we investigate how the current limits 
on CP violation in the decay D iiT^Tr+Tr" affect the 
model-independent measurement of 7 using the current 
approach in which CP conservation is assumed in charm 
decays. 

The first study placing limits on direct CP violation in 
D — > KgTT^Tr~ decay has been performed by the CLEO 
collaboration j 22| and has recently been improved by 
CDF [23]. Both measurements use isobar formalism to 
parametrize the D — > i^T^Tr+Tr" decay amplitude and as- 
sume that a CP asymmetry may appear in any of the 
quasi two-body amplitudes. Specifically, the amplitude 
is represented as the sum of resonance components for 
both the and D'^ decays: 

Ad =aoe'*° + ^ A+Mj , 

_ — (7) 

Ad ^aoe^'" +J2^J^i^ 
3 

where Ad and Ad are the amplitudes of and 
decays, respectively, hq and 5q are the amplitude and 
the phase of the non-resonant component (assumed to 
be CP-conserving), M.j and M.j are the quasi two-body 
resonant matrix elements (typically the relativistic Breit- 
Wigner amplitudes), and A"^ are the (complex) ampli- 
tudes of the resonant components. In the CP-violating 
case, A^j ^ Aj . In general, the presence of a CP- 
violating amplitude can result in the difference of both 
the magnitudes and the phases of the complex numbers 
A'j and Aj . The parametrization adopted by CLEO and 
CDF is 

4=a,e'(^^-±*^)(l±6,/a,), (8) 

where aj and Sj are CP-averaged parameters, while bj/aj 
and (j^j are small parameters related to the CP violation. 

We perform Monte Carlo (MC) studies to estimate 
the systematic uncertainty in the 7 measurement aris- 
ing from the current limits on CP violation in the D 
KgTT^TT'' decay obtained by CDF. A parametrization of 



FIG. 1. Bias of 7 measurement in the fit without accounting 
for CP asymmetry, using the amplitude generated without 
CP asymmetry (filled square), as well as with CP asymmetry 
in the magnitude b/a = 10% (filled circles), and in phase (p = 
0.1 (open circles) in each ampUtude component. Statistical 
errors are comparable with the size of markers. 



the form ([5]) is used. The D KgTT^n^ amplitude 
model used to generate event samples is based on the 
Belle measurement [2jj. A CP asymmetry of 10% is in- 
troduced one-by-one for the amplitude and phase in each 
partial amplitude (i.e. bj/aj = 0.1 and 0j = 0.1 radian). 
Large number of flavor-tagged D, Dq'p, correlated DD, 
and DK^ samples are generated so that the sta- 

tistical error of 7 measurement does not exceed 0.2°. The 
values 7 = 70°, tb = 0-1, Sb — 130° are used at the gen- 
eration stage. The samples are then fitted to extract 7 
value and other related parameters without taking CP 
violation into account. The resulting values of 7 for each 
variation of the D'^ decay amplitude are shown in Fig. [T] 

After the bias in 7 for a 10% CP asymmetry is ob- 
tained from MC, we recalculate the bias associated with 
current experimental limits by taking the central values 
of CP asymmetries and their errors from the CDF mea- 
surement [2^ assuming a linear dependence of 7 bias on 
CP- violating parameters bj/aj and 4>j. The resulting 
contributions from CP violation in each resonance are 
given in Table HI Finally, we calculate the total error by 
summing up the central values of bias linearly for each 
contribution, and the errors quadratically. The errors 
of the individual components are assumed to be uncor- 
related. As expected, the resulting 7 bias is consistent 
with zero within its error since no evidence of CP viola- 
tion has been found by CDF. Thus the 7 uncertainty is 
taken as the error of the bias and amounts to around 3° . 
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TABLE I. Contributions of CP violating amplitudes in D — > 
Kg-K^Tv" decay measured by CDF 23] to the 7 measurement 
bias for each contributing resonance, and the total 7 bias. 



Resonance 


Contribution to 7 bias (°) 


Amplitude 


Phase 


CF a:* (892) 


+0.09 ± 0.27 


-0.87 ± 2.09 


CF K^iUSO) 


-0.05 ± 0.05 


-0.23 ± 0.35 


CF Jf2*(1430) 


+0 07 ± 12 


— 04 ± 07 


CF if* (1410) 


+0 01 + 02 


—0 21 ±0 37 




+0 27 ± 89 


— 24 ± 97 




—0.32 ± 0.21 


-0.25 ± 0.36 


fo(980) 


—0 02 + 1 3 


— 02 + 38 


f9fl270) 


-0.09 ±0.10 


-0.06 ± 0.09 


/o(1370) 


-0.09 ± 1.06 


±0.01 ±0.26 


p(1450) 


-0.02 ±0.19 


-0.09 ±0.22 




-0.31 ±0.78 


-0.09 ±0.62 


0-2 


-0.07 ±0.08 


-0.04 ±0.56 


DCS if* (892) 


-0.04 ±0.24 


±0.22 ±0.15 


DCS 7^0* (1430) 


±0.23 ± 0.44 


-0.12 ±0.21 


DCS 7^2* (1430) 


-0.30 ± 0.56 


±0.03 ± 0.04 


Total 


-2.65 ±3.17 
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FIG. 2. Bias of 7 measurement in the fit with CP asymmetry 
accounted for, using the amplitude generated without CP 
asymmetry (filled square), as well as with CP- violation in the 
magnitude b/a — 10% (filled circles), and in phase </!> — 0.1 
(open circles) in each amplitude component. Note different 
vertical scale compared to Fig. [T] 



IV. MODEL-INDEPENDENT ANALYSIS WITH 
CP VIOLATION IN D ^ K'^tt+tt- 

Here we present results of the MC study performed to 
estimate how the fit procedure which allows for CP vio- 
lation in Z3 — >■ Kg7T~^TT~ decay described in Section Ullaf- 
fects the statistical precision of the 7 measurement com- 
pared to the CP-conserving case. 

We have performed MC simulation with 10^ events of 
flavor-specific D — > KgTT^n^' of each flavor, 120000 cor- 
related ijj{3770) DD decays, 120000 Dcp ^ A'^tt+tt" 
decays, and 60000 DK^, D -J> K%tt+tt- decays 

of each B sign. This B sample size corresponds roughly 
to the data sample expected in the upgraded phase of the 
LHCb experiment and at the Super B factory. The ratio 
of ?A(3770) and B samples was taken to be the same as in 
current analyses; we expect that the sufficient sample of 
^'(3770) decays will be collected by BES-III experiment 
and future tau-charm factory. In addition to this sample 
denoted by the factor fc = 1, we repeat the simulation 
with four times smaller (fc — 1/4) and four times larger 
{k = 4) samples to check how the error scales with the 
sample size. Current world-average values for the param- 
eters of i? — ?> DK decays are taken: 7 = 70°, = 0.1, 
5b = 130°. A total of 1000 pseudoexperiments are gen- 
erated and fitted for the sample of each size. 

Each MC sample is fitted with two techniques: a) the 
one which assumes CP conservation \n D K^tt^tt^ 
decay, and b) the one which allows for CP violation, as 
described in Section|lTl Using the similar variations of the 
D — >■ ii'gTr+Tr" amplitude involving 10% CP asymmetry 
in each resonance component as in Section Hill we show 



TABLE II. Comparison of the precision of the model- 
independent 7 measurement for the fit procedures with and 
without accounting for the CP violation in D — > A'^tt^tt" 
decay. The size of the simulated sample defined by the factor 
k is described in the text. 



k 


a(7) (No CP, °) 


cr(7) (with CP, °) 


Ratio 


1/4 


2.932 ±0.081 


3.021 ±0.084 


1.030 ± 0.040 


1 


1.525 ±0.042 


1.612 ± 0.049 


1.057 ± 0.043 


4 


0.713 ±0.019 


0.775 ± 0.019 


1.088 ± 0.039 



that the bias in 7 measurement is consistent with zero 
(see Fig. [2]). We also compare the statistical precision 
of the two approaches from the spread of fitted 7 values 
between pseudoexperiments. As shown in Table [ill the 
reduction of statistical precision due to increased number 
of free parameters does not exceed 10%. 



V. CONCLUSION 

We have shown that the current best limits on CP 
violation in D ^ X^tt+tt" decay coming from the mea- 
surement performed by CDF [23| translates to system- 
atic uncertainty in the determination of the CKM phase 
7 from B^ — >• DK^, D K^tt+tt" decay of the order 
of 3°. While the current world-average precision of 7 is 
9 — 12° [2^, [2^ and is not limited yet by this uncertainty, 
the data sample to be collected by LHCb experiment be- 
fore its upgrade should allow measurement with a preci- 
sion around 5° in which the B^ — s- DK^, D — >■ KgTT~^TT~ 
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Dalitz analysis will have significant weight [27| . It is thus 
useful to study the CP asymmetry in Z? — >■ ifg7r+7r~ with 
a larger data sample {e.g. at B factories and LHCb) to 
reduce this uncertainty. 

In addition, we have shown that even if the D 
KgTT~^TT~ decay is found to exhibit CP violation, it is 
possible to account for it and perform the unbiased mea- 
surement of 7 in DK^, D — > Kgn^n^ decay 
in a model-independent way. Compared to the model- 
independent technique which assumes CP conservation 
in _D ^> Kg7r~^TT~ [3 [2l|, this method has more free pa- 
rameters which, however, leads to a reduction of the sta- 



tistical precision not exceeding 10%. 
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